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hoop  stress  considerably  which  has  an  adverse  effect  on  bore  crack 
propagation.  This  study  considers  the  stress  intensity  factors  due  to  a 
radial  crack  taking  the  Bauschinger  effect  into  consideration. 

The  elastic-plastic  interfaces  during  loading  and  unloading  in  the 
autofrettage  process  divide  the  tube  into  three-ring  regions.  The  residual 
stress  distribution  in  each  region  is  quite  different.  When  a  crack  grows 
from  one  region  into  another,  the  previous  method  using  functional  stress 
intensity  fails.  A  new  method  is  used  to  obtain  stress  intensity  factors  for 
a  radial  crack  growing  out  of  the  reverse  yielding  zone.  This  approach  is 
based  on  crack  face  weight  functions  obtained  by  Sha  using  stiffness 
derivative  finite  element  techniques  coupled  with  singular  crack-tip  elements. 
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INTRODUCTION 


An  early  brittle  failure  of  a  cannon  tube  during  the  1960's  prompted  a 
renewed  interest  in  fracture  analysis  of  cannon  pressure  vessels  (ref  1). 

After  careful  investigation  of  the  cause  of  fracture,  some  basic  design 
changes  were  made  to  prevent  any  further  failures.  One  of  the  design  changes 
was  to  introduce  a  compressive  residual  stress  near  the  inner  radius  of  the 
cannon  by  an  autofrettage  process.  Considerable  efforts  have  been  made  to 
predict  the  residual  stress  distribution  in  an  overstrained  tube  based  on 
various  material  models  (refs  2-5).  In  order  to  confirm  the  theoretical 
predictions  of  residual  stresses,  experimental  methods  (refs  6,7)  have  been 
initiated  in  our  laboratory. 

Since  the  current  cannon  tube  design  uses  a  wall  ratio  close  to  two,  the 
maximum  compressive  residual  stress  at  the  bore  for  a  fully  autof rettaged  tube 
is  about  85  percent  of  the  yield  stress  0o  of  the  material.  Most  of  the 
earlier  predictions  of  residual  stresses  were  based  on  the  assumption  of 
elastic  unloading.  According  to  Milligan  et  al  (ref  8),  the  high  strength  gun 
steel  has  a  very  high  Bauschinger  effect.  To  account  for  the  compressive 
yield  strength  reduction  due  to  the  Bauschinger  effect,  Underwood  and  Kendall 
(ref  9)  have  tried  to  estimate  the  residual  stress  distribution  and  its  effect 
on  fatigue  crack  growth  rate.  Recently,  Chen  obtained  a  closed  form  solution 
of  residual  stresses  in  autof rettaged  tubes  based  on  a  theoretical  material 
model  taking  both  Bauschinger  and  hardening  effects  into  consideration  (ref 
10).  The  new  residual  stress  distribution  is  very  different  from  that 


References  are  listed  at  the  end  of  this  report. 


obtained  earlier  based  on  the  assumption  of  elastic  unloading.  The  functional 
stress  intensity  method  developed  for  the  computation  of  stress  intensity 
factors  of  a  radial  crack  initiating  from  the  inner  radius  of  a  tube  may  fail 
due  to  the  presence  of  a  reverse  yielding  region  near  the  bore.  The  objective 
of  this  report  is  to  develop  a  numerical  method  to  overcome  such  a  difficulty 
and  to  study  the  Bauschinger  effect  on  stress  intensity  factors  for  a  radially 
cracked,  partially  autof rettaged  gun  tube. 

THE  BAUSCHINGER  EFFECT  AND  RESIDUAL  STRESSES 

The  phenomenon  that  a  material  lowers  its  elastic  limit  in  compression 
(tension)  subsequent  to  a  previous  stressing  in  tension  (compression)  beyond 
the  elastic  limit  is  called  the  Bauschinger  effect.  A  quantity  representing 
the  magnitude  of  the  Bauschinger  effect  is  the  Bauschinger  effect  factor 
(BEF).  This  is  defined  as  the  ratio  of  the  yield  stress  upon  reverse  loading 
to  the  initial  yield  stress  (oQ).  The  BEF  (f)  is  a  function  of  percent  over¬ 
strain  (GP).  The  graph  of  BEF  versus  percent  overstrain  obtained  by  Milligan 
et  al  (ref  8)  for  a  modified  4330  steel  having  a  martensitic  structure  is 
shown  in  Figure  1  which  was  used  by  Chen  in  his  computations  of  residual 
stresses. 

Taking  the  Bauschinger  effect  factor  f  into  consideration,  Chen  obtained 
the  closed  form  solution  of  residual  stresses  in  autof rettaged  tubes  (ref  10). 
He  assumed  a  material  model  which  exhibits  the  stress-strain  curve  shown  in 
Figure  2  during  tensile  loading  and  unloading  after  overstrain.  The 
assumption  of  elastic-perf ectly  plastic  loading  was  supported  by  the  fact  that 
very  little  strain-hardening  was  observed  in  the  tensile  test.  A  bilinear 


2 


model  for  elastic-plastic  unloading  was  assumed  since  a  large  slope  of 
strain-hardening  (m'E)  did  develop  after  the  occurrence  of  reverse  yielding. 


Referring  to  the  point  0*  in  Figure  2  as  the  origin  of  a  new  (primed) 
coordinate  system  (a',e'),  the  reverse  yielding  curve  can  be  expressed  as 

o' 


1  +  f  +  m' C/(l-m’) 


(1) 


where  C  *  (E/o0)e'?,  E  is  Young's  modulus,  and  e*  is  the  plastic  strain  in 
the  primed  coordinates.  The  final  residual  stress  state  (denoted  by  a  double 
prime)  is  obtained  by  summing  the  stress  state  corresponding  to  the 
elastic-perf ectly  plastic  loading  and  the  primed  stress  state  corresponding  to 
unloading.  The  tangential  and  radial  components  can  be  written  as 

o9"  =  o9  +  00’  ,  or"  =>  or  +  or'  (2) 

Elastic-Plastic  Loading 

Let  a  and  b  be  the  inner  and  outer  radii  of  the  cylinder,  respectively, 
and  let  the  material  be  elastic-perf ectly  plastic,  obeying  Treses’ s  yield 
criterion.  The  stress  components  are  given  in  Reference  2  as  follows: 
ar/0o 

°0/oo 


1  p  P 

-(+l  +  --)-to-  ,  a  <  r  <  p 

2  b2  r 


(3) 


1  p£  P2 
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ar/°o 

a8/°o  2  b‘ 

where  P  is  the  elastic-plastic  boundary  relating  to  the  internal  pressure  p  by 


p  1  Pz  P 

—  =»-(!-  --)  +  In  (-) 


(3) 


The  equivalent  plastic  strain  can  be  calculated  by 


E  Pz 

“  ep  *  8l  (-o  *  1)  ,  a  <  r  <  p 


where 


h  -  —  (i-\>2) 
✓5 


In  Eq.  (7)  v  is  Poisson's  ratio. 


Elastic-Plastic  Unloading 


If  the  internal  pressure  p  is  subsequently  removed  completely,  the 


unloading  may  be  either  elastic  or  elastic-plastic  depending  on  whether  the 


magnitude  of  p  is  less  than  or  greater  than  p^  which  is  the  minimum  pressure 


to  cause  the  reverse  yielding  to  occur.  If  p  <  pm,  the  unloading  is  entirely 


elastic  and  the  stress  components  are 


°r  p  b^ 

- - [±  --  -  1] 

00’  b2  r 

~  -  1 


Using  Tresca's  criterion  subject  to  cr"  >  oz"  >  o0”,  the  reverse  yielding  will 


not  occur  if 


or"  -  00"  <  fo0 


Substituting  Eqs.  (2),  (3),  and  (8)  into  Eq.  (9),  we  obtain  the  expression  for 


Pm  1  a4 


-  -  (l+f)(l  -  --) 


ao  2 


If  p  >  pm,  reverse  yielding  will  occur  in  a  region  a  <  r  <  p'  with  p'  <  p 


upon  unloading  and  we  have  from  Eq.  (1) 


V  -  OQ"  ,  fo0  +  m’Ee'p/d-^)  ln  a  <  r  <  p' 


The  stresses  in  the  plastic  and  elastic  regions  are  given  by  Bland  (ref  3)  as 

1  p'  a2  r 

V/o0  -  P /o0  -  -  02'a+f)(")2(l  ~  -5)  -  (l-B2')(1+n«n(-)  ,  (a<r  <p' )  (12) 

z  ^  r  Si 


o9'/a0  =  or’/o0  -  (1+f)  -  m' C / (1-m' )  ,  (a<r<p')  (13) 

°r'/ao  1  p'  ,  p'  ~  (^) 

=  -  d+f)[i(")2  -  (— )2]  ,  ( P  <r  <b) 

o0'/ao  2  r  b  (15) 

where 

8l'  -  (l-m')/[m'+(l-m')/01]  (16) 

02*  -  m'  * / (l—m*  )  (17) 

C'  *  8i '  (l+f)(  p'  2/r^*l)  (18) 


The  value  of  p'  can  be  found  from  continuity  of  or'  at  p'  from  Eqs.  (12)  and 
(14).  Figure  3  shows  different  residual  hoop  stress  distributions  for  a 
hollow  cylinder  with  a  wall  ratio  of  two  and  E/aQ  ■  200,  v  ■  0.3,  p/a  “  1.6 
and  2.0.  The  broken  lines  are  elastic  unloading  solutions  (f«l),  while  the 
solid  lines  are  elastic-plastic  unloading  with  m*  *  0  and  0.3.  The  figure 
shows  a  drastic  difference  in  oq"  in  the  reverse  yielding  region  near  the  bore 
and  it  shows  that  p'  varies  slightly  with  P  and  m' . 

FUNCTIONAL  STRESS  INTENSITY  METHOD 

We  have  developed  the  functional  stress  intensity  method  (ref  11)  for  the 
computation  of  stress  intensity  factors  of  radial  cracks  in  a  pressurized  and 
partially  autof rettaged  cylinder  by  combining  the  finite  element  method  and 
the  weight  function  method.  A  weight  function  vector  h  is  a  universal 
function  which  depends  only  on  geometry  and  not  on  loadings  (ref  12).  For  a 
given  radially  cracked  ring,  if  Kj(l)^  the  mode  I  stress  intensity  factor,  and 


W-, 

'M 


vd),  the  normal  component  of  crack  face  displacement  associated  with  a 
symmetric  load  system  1,  are  known,  then  the  normal  component  of  the  crack 


face  weight  functions  can  be  expressed  by 


hly(x > 


2 


where  H  *  E  for  plane  stress  and  H  “  E/(l-v2)  for  plane  strain,  x  =*  (r-a)  is 
distance  measured  from  the  base  of  the  crack  along  the  crack  face  toward  the 
crack  tip,  and  l  is  the  crack  length.  For  any  symmetric  load  applied  to  the 
same  cracked  ring,  the  stress  intensity  factor  associated  with  the  new  load 
can  be  found  from 


1 


H  l  3v(L) 

---  J  p(x)  -----  dx 


Kl(l)  o  3 l 

where  p(x)  is  the  normal  stress  in  the  tangential  direction  at  the  crack  site 
due  to  the  new  load  applied  to  the  uncracked  ring.  If  the  uniform  tension  p0 
at  the  outer  radius  is  taken  as  load  one,  we  can  find  numerical  values  of 
and  by  the  finite  element  method.  We  wish  to  use  Eq.  (20)  to 

obtain  Kj  associated  with  the  autofrettage  residual  stresses.  It  can  be  seen 
from  Eqs.  (2),  (3),  (8),  and  (12)  through  (15)  that  the  residual  stress  Oq" 
has  the  general  expression 

°9”(r)  =»  A]_  +  ki/r2  +  A3in(r)  (21) 

for  a  tube  subject  to  an  elastic-plastic  loading  prior  to  an  elastic 
unloading,  where  A^,  i  -  1,2,3  are  superposition  constants. 

Let  p(x)  in  Eq.  (20)  be  Oq”  above.  The  integration  to  find  Kj  requires 
an  expression  for  3v(l)/3£,,  A  method  which  assumes  v(l)  as  a  conic  section 
(ref  13)  has  been  used  by  Grandt  (ref  14).  Another  method  which  avoids  the 
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need  of  3v(*)/3£  has  been  developed  by  Pu  (ref  11).  He  uses  finite  elements 
to  obtain  functional  stress  Intensities  Kc(l),  ^(r-2),  and  Kc(  ter)  defined  by 


H  X  3v(^) 

Kc(p)  =* - /  p - dx  (22) 

C  K(l)  o  3A 

for  a  crack  face  loading  p  *  1,  p  *  r“2,  and  p  ■  te(r),  respectively.  Once 
the  functional  stress  intensities  are  known,  then  the  stress  intensity  factor 
Kj(p)  can  be  found  for  a  general  residual  stress  distribution  of  Eq.  (21)  by 
an  algebraic  equation 

Kr(p)  -  A1KC(1)  +  A2Kc(r-2)  +  A3KC( ter)  (23) 

This  method  was  successfully  applied  to  various  residual  stresses  predicted 
from  various  material  models  (ref  15). 

Nc-  for  a  radial  crack  in  a  reverse  yielding  zone,  the  stress  intensity 
factor  can  still  be  found  by  the  functional  stress  intensity  method  since  the 
residual  stress  remains  the  same  form  as  Eq.  (21)  with  the  superposition 
constants  obtainable  from  Eqs.  (3)  and  (13).  However,  if  the  radial  crack  is 
longer  than  the  reverse  yielding  zone,  the  functional  stress  intensity  method 
fails  since  the  crack  face  loading  p(x)  has  two  different  expressions  in  two 
different  regions:  0  <  x  <  p'-a  and  p'-a  <  x  <  £.  Let  py(x)  and  pe(x)  be  the 
crack  face  loading  (residual  hoop  stress)  in  the  reverse  yielding  region  and 
the  elastic  region.  Equation  (20)  becomes 

,p’-a  X 

KI  “  2 /  PY(x)hIy(x, *)dx  +  2 /  Pe(x)hIy(x,  £)dx  (24) 

o  3  p'-a  3 

To  use  this  equation,  we  have  to  find  the  explicit  crack  face  weight  function 
hjy  using  the  stiffness  derivative  finite  element  technique  explained  in  the 


next  section. 


EXPLICIT  CRACK  FACE  WEIGHT  FUNCTION 

The  stiffness  derivative  method  for  determining  linear  elastic  stress 
intensity  factors  was  introduced  by  Parks  (ref  16) .  The  method  was  extended 
to  calculate  the  weight  function  vector  field  by  Parks  and  Kamenetzky  (ref 
17).  An  efficient  finite  element  evaluation  of  explicit  weight  functions  has 
been  established  by  Sha  (ref  18)  who  combined  the  stiffness  derivative 
technique  with  special  singular  crack-tip  elements.  He  and  Yang  applied  the 
method  to  radial  crack  problems  of  a  hollow  disk  (ref  19).  The  degenerated 
quarter-point  quadratic  elements  were  used  around  a  crack  tip.  These  singular 
elements  were  surrounded  by  the  standard  eight-node  quadrilateral  elements. 

The  virtual  crack  extension  of  an  amount  6i  was  simulated  by  advancing  the 
crack-tip  node  by  in  the  direction  colinear  with  the  mode  I  radial  crack 
(x-direction) .  The  surrounding  quarter-point  nodes  were  also  shifted  to  new 
locations  and  there  was  no  change  in  location  of  all  other  nodes.  Hence  only 
a  few  crack-tip  elements  have  experienced  changes  in  elemental  stiffness  due 
to  the  virtual  crack  extension.  This  makes  the  stiffness  derivative  technique 
very  efficient  from  a  computational  viewpoint. 

From  the  displacement  form  of  the  finite  element  method,  the  equation  of 
equilibrium  is 

[K]  {u}  -  {F}  (25) 

where  [K]  is  the  global  sitffness  matrix,  and  (u)  and  {F}  are  displacement 
vector  and  load  vector,  respectively.  The  change  in  displacement  per  unit 
crack  extensions  can  be  found  by  differentiating  Eq.  (25)  with  respect  to 


crack  length  i 


Mote  chat  d{F}/d£  *  0  if  we  select  a  load  system  which  consists  of  only 
surface  tractions  not  applied  on  the  crack  face.  The  global  stiffness 
derivative  d(K]/di  is  the  sum  of  Nc  elemental  stiffness  derivatives 


d[K]  Nc  ^[kiJ 
d  Z  d  i 


(27) 


where  Nc  is  the  number  of  crack-tip  elements  and  [kjJ  is  the  element  stiffness 
of  the  ith  crack-tip  element.  For  a  small  crack  extension  61,  the  element 
stiffness  derivative  may  be  approximated  by 


d[ k.±]  [kil  £+6£  "  tkil  % 
dl  61 


(28) 


where  [kjj£+$£  and  [k^Jj,  are  element  stiffness  after  and  before  virtual  crack 
extension,  respectively. 


The  displacement  vector  is  a  function  of  position  (x,y)  and  the  crack 
length  1  as  follows: 

{u}  -  (u(x,y,l)}  (29) 

Applying  the  chain  rule  of  differentiation  of  Eq.  (29)  with  respect  to  1 
gives 


3  (u  }  d  (u  }  3  {u  }  dx  3  (u  }  dy 

—  — — —  3  — — — —  —  mmmm  mm  m  MMMM  mm 

31  dl  3x  dl  3y  dl 


(30) 


For  a  mode  I  crack  lying  on  the  x-axis  with  colinear  virtual  crack  extension, 
we  have  dx/dl  ■  1  and  dy/d 1  •  0.  Definitions  of  global  coordinates  and 
displacements  for  isoparametric  elements  are 


u 


(31) 


x  -  I  ,  u  -  l  Ni(C,n)ui 

y  m  l  Ni(5,n)yi  ,  v  -  l 
where  5,h  are  local  coordinates;  x,y  are  global  coordinates;  x^,y^  are  global 
coordinates  of  node  i;  and  are  x  and  y  components  of  displacement  at 
node  i;  N^(5,h)  are  shape  functions  which  interpolate  the  displacement  over 
the  element.  The  finite  element  evaluation  of  3{u}/3x  can  be  carried  out  from 

3{u>  l  £y  3tNiuil  3y}  32 

3x  det[J]  35  3h  3n  3? 

where  [J]  is  the  Jacobian  matrix.  This  leads  to  the  expression  for  the  y 
component  of  mode  I  weight  function  vector  at  (x,y) 

H  dv(!>  1  3[NiUi]  ay  ^1%^]  ^ 

Kj(^-)(£)  ^  di.  det[J]  3^  3n  35  35 

For  a  hollow  disk  of  b/a  *  2,  Sha  and  Yang  [19]  have  obtained  explicit 
weight  functions  for  a  single  bore  or  rim  radial  crack.  In  a  private 
communication,  Sha  has  provided  the  following  expression  to  approximate  the 
dimensionless  crack  face  weight  function  component 


2/S  hIy(~)  -  l  Dn(~) 
1  n-1  * 


n/2-i 


(34) 


where  r3,  the  distance  from  the  crack  tip  along  the  crack  face,  is  related  to 
x  and  r  by 

r  g  -  l  -  x  m  Jl+a-r  (35) 

The  coefficients  Dn,  n  ■  1,2, 3, 4,  determined  by  the  least  squares  technique, 
are  given  in  the  following  table  for  various  crack  lengths. 
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TABLE  1.  LEAST  SQUARES  FITTED  COEFFICIENTS  OF  EQUATION  (34)  FOR  A  SINGLE  BORE 


CRACK  IN  A  HOLLOW  CYLINDER  OF  b/a  -  2  WITH  VARIOUS  CRACK  LENGTHS  c 
*/(b-a) 


l 

b-a 

T 

1 

1 

I>1 

T 

1 

1 

d2 

T 

1 

1 

1 

°3  i 

d4 

0.01 

T 

1 

I 

0.79660+00 

T 

1 

1 

-0.77210-02 

T 

1 

| 

0.29350+00  | 

0.36200+00 

0.02 

1 

1 

i 

0.7976D+00 

1 

1 

i 

-0.6980D-02 

1 

1 

i 

0.2840D+00  1 

0.3608D+00 

0.03 

1 

1 

1 

0.8130D+00 

1 

1 

| 

-0.7176D-01 

1 

1 

1 

0.37110+00  | 

0.3134D+00 

0.05 

i 

1 

i 

0.8008D+00 

1 

1 

i 

-0. 5811D-02 

i 

1 

| 

0.26490+00  | 

0.3551D+00 

0.06 

1 

i 

0.83350+00 

1 

1 

i 

-0.17140+00 

I 

1 

i 

0.49420+00  | 

0.24210+00 

0.08 

1 

1 

1 

0.80510+00 

1 

1 

i 

-0.2789D-01 

1 

1 

0.31160+00  | 

0. 3270D+00 

0.10 

1 

1 

i 

0.8286D+00 

1 

1 

i 

-0.1818D+00 

1 

1 

i 

0.61190+00  | 

0.1554D+00 

0.20 

1 

1 

i 

0.81320+00 

1 

1 

i 

-0.8254D-01 

1 

1 

i 

0.57690+00  | 

0.24140+00 

0.30 

1 

i 

0.80260+00 

1 

1 

i 

0.1409D-02 

1 

1 

i 

0.58810+00  | 

0.37560+00 

0.40 

1 

1 

i 

0.79810+00 

1 

1 

i 

0.4661D-01 

1 

1 

| 

0.6908D+00  | 

0.49350+00 

0.50 

1 

i 

0.79560+00 

1 

1 

i 

0.7236D-01 

1 

1 

i 

0.9065D+00  | 

0.5437D+00 

0.60 

1 

i 

0.79900+00 

1 

i 

0.4395D-01 

1 

1 

i 

0. 1361D+01  1 

0.4059D+00 

0.70 

1 

1 

i 

0.80040+00 

1 

1 

i 

0.2255D-01 

1 

1 

i 

0.21300+01  | 

-0.3915D-01 

0.80 

1 

1 

JL 

0.79760+00 

1 

1 

_L 

0.4388D-01 

1 

1 

_L 

0.3537D+01  I 

-0.11310+01 

STRESS  INTENSITY  FACTORS 


To  use  Eqs.  (34)  and  (24),  we  need  a  certain  transformation  of  variables. 
Let  us  use  a,  the  bore  radius,  to  normalize  all  linear  lengths  and  use  the 
same  notations  (before  normalization)  to  denote  the  normalized  lengths  except 
that  a  is  unit  and  b  is  w,  the  wall  ratio.  Denoting  ra/ l  by  t  in  £q.  (34), 
then  Eq.  (24)  associated  with  the  residual  hoop  stress  becomes 


Kr  i  Oq"  4  t*  Og"  4 

- -  -  /  ( - )Y  (  I  DnW2-1)dT+  /  ( - )e(  l  DnTn/2-l)dT 

0-/I  T'  °o  n«l  0  °o  n*l 


1  Oq"  4  i  Oq"  Oq"  4 

-  /  ( - >e  (  l  DnW2-l)dT  4-  /\  [( - )  -  ( - )e]  [  £  DnW2-l]dT  (36) 

°  °o  n“l  T  °o  Y  0o  n-1 

where  (0g"/oo)e  is  the  residual  hoop  stress  in  the  elastic  unloading  region, 

while  (Oq"/o0)y  is  that  in  the  reverse  yielding  region  and  t'  is  the  value  of 

x  corresponding  to  r  *  p',  the  elastic-plastic  interface  during  unloading. 

The  residual  hoop  stress  can  be  represented  by  the  general  form 

°0" 

( - )e  *  Aie  +  A2e  r_2(t)  +  A3e  *n(r/x)  (37) 

ao 


°0"  , 

(  )y  a  Aiy  +  A2y  r”2(x)  +  A3Y  £n(r/x)  (38) 

oQ 


where  the  coefficients  Aie . A3Y,  can  be  obtained  from  Eqs.  (2),  (3),  and 

(12)  through  (15).  From  Eq.  (35)  and  the  definition  of  T,  we  have 

r(T)  -  *(  VT)  .  *o  “  (1+*)/*  (39) 

Denoting  Iin(a,8),  l2n(a»0)»  a^d  ^3nC°»0)  as  the  following  definite  integrals, 
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Iln(a,S)  -  /0  W2“ldr 
a 

l2n («.e)  -  /  r-2(T)W2-ldT  a  =  1,2, 3, 4  (40) 

a 

l3n(“,e)  "  J*  fa(r(t))W2-ldT 
a 

Eq.  (36)  can  be  written  in  the  fora 

Kr  i  3  4  3  4 

Aie  I  Vi^.D  +  I  (AiY-Aie>  I  DnIin(  f ,  1)  (41) 

a0/nT  /it  i=«l  n*l  i”l  n*l 

The  expressions  of  integrals  I^n(  o,  S) ,  i  *  1,2,3,  n  *  1,...4,  are  given 
in  the  Appendix.  Equation  (41)  can  be  used  for  any  loading  which  yields  a 
hoop  stress  of  the  fora  given  by  Eq.  (21)  and  for  any  multiple-cracked 
cylinder  as  long  as  the  coefficients  Dn  are  known  for  that  particular  cracked 
geometry.  In  case  there  is  no  reverse  yielding  region,  Eq.  (41)  reduces  to 

KI  1  3  4 

-  “  ”  "  I  Ale  l  DnIln(0,l)  (42) 

Og/ttfc  Vt[  i»l  n**l 

For  a  shallow  crack  which  lies  entirely  in  the  reverse  yielding  region,  Eq. 
(41)  becomes 

KX  ]_  3  4 

l  AiY  l  DnIin(0,D  (43) 

o0/tt l  /it  i=l  nal 

Table  1  gives  Dn  for  a  single  bore  crack  of  various  discrete  crack  lengths. 

For  a  crack  length  not  given  in  Table  1,  Sha  and  Yang  (ref  19)  suggested  the 
use  of  cubic  spline  interpolation  to  obtain  the  weight  function  for  that 
length  from  discrete  values  in  Table  1  followed  by  the  least  squares  technique 
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Co  calculate  D's. 


Equation  (42)  is  used  to  compute  Kj/p0/iT£  for  a  single  bore  crack  in  a 
hollow  cylinder  of  w  *  2  subject  to  uniform  tension  pQ  at  the  outside 

cylindrical  surface.  The  Lame'  solution  gives  »  A2  *  w2/(w2-l)  and  A3  »  0. 

Numerical  results  for  various  crack  lengths  are  very  accurate  in 
comparison  with  results  previously  reported  in  Reference  20.  A  further  check 
of  numerical  computations  of  Kj/a0/itT  from  Eq.  (42)  has  been  done  for  a 
partially  autof rettaged  cylinder  without  reverse  yielding.  The  results 
confirm  those  published  in  Reference  21. 

For  a  cylinder  of  w  *  2,  the  Bauschinger  effect  factors  for  60  percent, 

80  percent,  and  100  percent  overstrain  are,  from  Figure  1,  f  ■  0.44,  0.42,  and 
0.38,  respectively.  Using  v  ■  0.3,  E/ aQ  »  200,  the  elastic-plastic  interface 
p'  during  unloading  depends  on  f  and  m' .  It  varies  from  p'  *  1.106  for  e  ■ 

0.6  (f  *  0.44)  and  m'  *  0.3  to  p'  ■  1.20  for  e  -  1  (f  ■  0.38)  and  m'  ■  0.  The 

residual  hoop  stress  distributions  near  the  bore  for  these  two  cases  are  shown 
in  solid  lines  in  Figure  3  while  broken  lines  are  corresponding  stresses 
without  taking  Bauschinger  effect  into  consideration.  The  residual  stress  in 
the  reverse  yielding  region  varies  drastically  with  f  and  m' .  Stress 
intensity  factors  due  to  residual  stresses  are  calculated  from  Eq.  (43)  or 
(41)  depending  on  whether  the  crack  tip  is  inside  or  outside  of  the  reverse 
yielding  zone.  Figures  4  and  5  are  graphs  of  Kj/(o0/irl)  as  a  function  of 
dimensionless  crack  length  c  ■  £/(w-l)  for  e  *  0.6  and  e  ■  1,  respectively. 
Superposing  stress  intensity  factors  due  to  an  internal  pressure  p  ■  o0/Lf, 
the  combined  stress  intensity  factors  from  both  p  and  residual  stresses  are 
shown  in  rigures  6  and  7  for  Lf  ■  3  and  for  £  •  0.6  and  e  ■  1.0,  respectively. 


Figure  8  is  a  similar  graph  for  Lf  ■  1.5  and  e  »  1.0.  In  Figures  6,  7,  and  8, 
the  curves  in  the  centerline  correspond  to  m'  »  0,  the  curves  in  the  broken 
line  are  for  m'  =  0.3,  while  the  solid  lines  are  results  obtained  earlier 
without  taking  reverse  yielding  into  consideration.  For  c  =  0.01,  Figure  7, 
values  of  Kj/(o0/irT)  are  0.55  and  0.37  for  m'  -  0  and  0.3,  respectively, 
versus  the  corresponding  value  0.058  from  the  solid  line.  This  indicates 
the  Bauschinger  effect  will  greatly  reduce  the  advantageous  effect  of 
compressive  residual  hoop  stress  introduced  by  the  autofrettage  process. 

CONCLUSION 

The  Bauschinger  effect  of  the  high  strength  gun  steel  causes  the  presence 
of  the  reverse  yielding  in  a  partially  autof rettaged  cannon  tube  of  wall  ratio 
of  two.  This  reverse  yielding  reduces  the  magnitude  of  compressive  hoop 
stress  considerably  in  the  reverse  yielding  region.  This  stress  reduction 
will  result  in  much  higher  stress  intensity  factors  in  a  shallpw  bore  crack  in 
a  pressurized  and  autof rettaged  gun  tube.  The  higher  the  stress  intensity 
factor  implies  the  lower  the  fatigue  life. 

The  expression  in  Eq.  (34)  for  explicit  crack  face  weight  function 
obtained  by  a  combination  of  stiffness  derivative  technique  and  special 
singular  crack-tip  elements  is  highly  accurate.  It  recovers  the  previous 
stress  intensity  factor  results  for  a  single  bore  crack  in  a  tube  subject  to 
various  loading  conditions.  The  method  can  treat  any  crack  face  loading 
including  stress  discontinuities  and  stress  gradient  discontinuities  over  the 


Numerical  results  here  are  limited  to  single  bore  radial  cracks. 
However,  the  method  is  general  for  either  bore  cracks  or  rim  cracks  for  any 
number  of  cracks  as  long  as  we  have  obtained  the  explicit  crack  face  weight 
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PERCENT  TENSILE  OVERSTRAIN 


Bauschinger  effect  factor  V9.  percent  tensile  overstrai 
martensitic  structure  of  a  4330  modified  steel. 


Figure  4.  Dimensionless  stress  Intensity  factors  as  a  function  of 

dimensionless  crack  length  c  due  to  compressive  hoop  stress  in  a 
cylinder  having  60  percent  degree  of  autofrettage  for  various 
values  of  f  and  m' . 
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Figure  5.  Dimensionless  stress  Intensity  factors  as  a  function  of 

dimensionless  crack  length  c  due  to  compressive  hoop  stress  In 
cylinder  having  100  percent  degree  of  autofrettage  for  various 
value 8  of  f  and  m' • 
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Figure  6.  Resultant  stress  intensity  factors  as  a  function  of  crack,  length  c 
for  various  values  of  f  and  o'  in  a  pressurized,  autofrettaged 
cylinder.  Internal  pressure  p  *  aQ/ 3  and  degree  of  autofrettage  e 
-  0.6. 
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Figure  7.  Resultant  stress  intensity  factors  as  a  function  of  crack  length  c 
for  various  values  of  f  and  o'  in  a  pressurized,  autof rettaged 
cylinder.  Internal  pressure  p  -  oQ/3  and  degree  of  autofrettage 
e  -  1.0. 
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Pigure  8.  Resultant  stress  Intensity  factors  as  a  function  of  crack  length  i 
for  various  values  of  f  and  n'  in  a  pressurized,  autof rettaged 
cylinder.  Internal  pressure  p  -  c0/i.5  and  degree  of  autofrettage 


APPENDIX 


The  definite  integrals  defined  by  Eq.  (40)  are  explicitly  given  below: 
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